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ABSTRACT 

We present a time-dependent cosmic-ray modified shock model for which the calculated Ha emis- 
sivity profile agrees well with the Ha flux increase ahead o f the Balmer-dominated shock at knot g in 
Tycho's supernova remnant, observed bv lLee et al.l (|2007l ). The backreaction of the cosmic ray com- 
ponent on the thermal component is treated in the two-fluid approximation, and we include thermal 
particle injection and energy transfer due to the acoustic instability in the precursor. The transient 
state of our model that describes the current state of the shock at knot g, occurs during the evo- 
lution from a thermal gas dominated shock to a smooth cosmic-ray dominated shock. Assuming a 
distance of 2.3 kpc to Tycho's remnant we obtain values for the cosmic ray diffusion coefficient, k, 
the injection parameter, e, and the time scale for the energy transfer, r, of k — 2 x 10 24 cm 2 s _1 , 
e = 4.2 x 10~ 3 , and r = 426 y, respectively. We have also studied the parameter space for fast 
(300 kms -1 < v s < 3000 kms^ 1 ), time-asymptotically steady shocks and have identified a branch of 
solutions, for which the temperature in the cosmic ray precursor typically reaches 2-6 x 10 4 K and the 
bulk acceleration of the flow through the precursor is less than lOkrns -1 . These solutions fall into 
the low cosmic ray acceleration efficiency regime and are relatively insensitive to shock parameters. 
This low cosmic ray acceleration efficiency branch of solutions may provide a natural explanation for 
the line broadening of the Ha narrow component observed in non-radiative shocks in many supernova 
remnants. 

Subject headings: acceleration of particles — hydrodynamics — line: profiles — methods: numerical 
— shock waves — supernova remnants 
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1. INTRODUCTION 

Balmer-dominated filaments in supernova remnants 
(SNRs) trace fast, non-radiative shocks that propagate 
into partially neutral, diffuse media. The filament s are 
sheets of shocked gas seen edge-on (jHesterl Il987f ) and 
have been observe d and studied in many SNRs includ- 



ing Tycho's SNR (fChevalier et all 119801: iKirshner et al 
19871 ISmith et all 1199 U iGhavamian etjaT "' 



2000, 



hereafter GOO and G01, 



_ , 2001 
respectively; I Lee et aU 12007 



hereafter L07), the Cyg nus Loop (e.g. iRavmon d et al 
1983]: iHester et al.lfl99l RCW86 (e.g. iSollerman et al 



20031 Kepler's SNR (e.g.lBjair_eVaL 1991; Sa nkrit et al 



20031 . SN1006 (e.g. lOhavamian etafl [200l. a nd four 
remnants in the LMC (e.g. iSmith et aU I1994D . This 
work focuses on modeling a particular Balmer-dominated 

shock located at knot g in Tycho's SNR (after 

iKamper fc van den Berghl [1978) and recently observed 
by L07. 

The Ha spectral line profile of a Balmer-dominated 
filament has narrow and broad components, whose 
widths represent the preshock temperature of neutral 
H and the postshock temperature of protons, respec- 
tively (jChevalier fc Raymond! Il978fl . The width of the 
broad component and the ratio of intensities are diag- 
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nostics for the shock speed and the degree of electron- 
ion temperature equilibration behind the shock. Shock 
models t o quantify thes e diag nostics were first devel- 
oped by iCheyalier et aU (1 19801) and later i mprov ed by 
Smith et aU (11991I L G01. iHeng fc McCravl (|2007h . and 
Currently, the most advanced non- 



Heng et al.l 



radiative shock models are thos e of |van Adelsberg et al.l 
(|2008f) . Except for the work bv lBoulares fc Coxl (|1988fT 
shock models used to interpret optical observations, to 
date, have not included modifications of the shock struc- 
ture due to diffusive shock acceleration (DSA) of cosmic 
rays (CRs). 

Balmer-dominated filaments are usually associated 
with forward shocks of the expanding SNR bub- 
bles. The only remnant in which B aimer emission 
from a reverse sho ck has been observed is SN 1987 
([Michael et al.ll2003l ). Forward shocks in SNRs are also 
sites of CR electron acceleration, as evidenced by syn- 
chrotron radio and X-ray emission (iKovama et al1 l!995 



Gotthelf et al.ll2001l:lLong et alj2003tlBamba et alj2005 



Cassam-Chenai et al.ll2007|). as well as likely sites of CR 
ion a c celeration (Blandford & Eichlcr 19871: iDrurv et ail 
120011: I Warren et all l2005t ). In most cases, Balmer- 
dominated shocks do not show evidence for synchrotron 
emission, which suggests that particle acceleration in 
Balmer-dominated shocks is not efficient. The neu- 
tral component in the upstream thermal gas that is 
required for non-radiative shocks to produce Balmer- 
dominated filaments may be damping the turbulence 
necessa ry for efficient cos mic ray acceleration in SNR 
shocks (jDrury et al.lll996l ). Conversely, the heating of 
the upstream medium due to efficient CR acceleration 
may prevent most neutrals reaching the shock before be- 
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ing ionized (|Hester et al.lll994j ). Two known exceptions 
where Balmer emission and synchrotron X-ray emission 
coincide are knot g in Tycho's rem nant and a small por- 
tion o f the eastern rim of SN 1006 (|Cassam-Chena'i et al.l 
2008). However, a direct connection between the Balmer 
emission producing shock and the X-ray producing shock 
cannot be made for either case because the X-ray mor- 
phology is not resolved to the level of the optical emis- 
sion and because the effects of projection are uncer- 
tain. Many Balmer-dominated filaments are observed to 
bound regions of X-ray emission whose spectra are con- 
sistent w ith thermal emission of a shock-heated am bient 
medium (jHester et al.lll994 iRavmond et al.ll2007| ). 

The theory of DSA predic ts a CR pr e cursor 
ahead of the gas su bshock (|Drurv fe Volkl 119811 ; 
iBerezhko fc Ell ison 1993), in which the upstream gas is 
pre-heated and accelerated over a characteristic distance 
k/v s , where k is a momentum averaged CR diffusion co- 
efficient, and v s is the shock speed. The value of k for CR 
ions depends on the spectrum of the magnetic wave-field, 
tho ught to be generated by the CR streaming instabil- 
ity (|Belll [P978). and has not been well constrained by 
observations yet. While in the general ISM n ~ 3-5 x 
10 28 cm 2 s _1 ([Strong et~aTll2007f) . at the forward shocks 
of SNRs k is thought to be close to the Bohm diffusion 
limit k b w 3 x 10 22 cm 2 s- 1 {fiBo/^Gy 1 (p/GeVc^ 1 ), 
where (3 is the ratio of particle speed to the speed of 
light, B is the large sca le magnetic field, and p is the 
CR particle momentum (|Drury|[l983f ). Kb is the lower 
limit to k allowed by the standard theory of DSA, and 
corresponds to saturat ed field fluctuations, SB /Bo = 1. 
ISollerman et al.1 (|2003l ) estimated upper limits for the dif- 
fusion coefficient of CR ions in several SNRs in the range 
k ~ 10 25 cm 2 s _1 -2 x 10 27 cm 2 s _1 from the condition 
that neutrals must survive ionization while experiencing 
the amount of heating implied by the Ha n arrow compo- 
nent l inewidth. The expression derived bv lParizot et al.l 
(2006) for the electron diffusion coefficient as a function 
of the synchrotron X-ray cutoff energy implies that the 
electron diffusion coefficient for the SNRs studied in their 
work is within only a factor of a few greater than the 
Bohm diffusion coefficient. 

The spectral shape and the sharply peaked radial pro- 
files of the synchrotron X-ray emission at the forward 
shocks of young SNRs require a post shock magnetic 
field s trength of the o r der of lOOuG (jVink fc Lamind 
[20031: iVolk et al.l [20051 : iBalletl I2006T) . Since compres- 
sion alone is insufficient to produce such a gain in field 
strength, it is thought that perturbations in the preshock 
medium generated by the CR streaming instability are 
non- l i nearly amplified beyond 5B/B n = 1 (jBell fc Lucekl 
120011: lBelll 12004 iParizot et al.l 12006). While some sim- 
ulations show that there exist rapidly growing nonres- 
onant modes (jZirakashvili et al.ll2008| ). t hese are satu- 
rated at 5B/Bq ~ 1 in other simulations (Niemicc et al. 
2008), and the problem of magnetic field amplification 
in the CR precursor remains an unresolved. Another 
mechanism for wave generation is the CR driven acous- 
tic in stability (iDrurvl 11984 iDrurv fc Falielll986t IChalovi 
[1985 iKang et all 11992). which may play a role in pro- 
longing the c onfinement of high energy CRs in the 
CR precursor (IBerez hko 1986 UMalkov fc Diamondll2006t 
iDiamond fc Malkovll2007ft . 



In general, strong MHD turbulence in the CR precur- 
sor leads to energy dissipation by wave damping that 
will affect the structure of th e shock (|MacKenzie fc Voelkl 
Il982t ICaprioh et a Tl l2008bl). Althou gh theories of wave 
damping exist (see e.g. IWhandll997l ). neither the rate of 
energy dissipation nor the fractions of the dissipated en- 
ergy going into internal energy of the various components 
of the flow (CR electrons and ions, and thermal elec- 
trons and ions) are known from observations. In models 
of DSA, the original treatme nt of Alfven wave damp- 
ing in CR modified shocks by IV51k fc McKei^ii (fl98lh 
is commonly adopted. The damping of acoustic waves, 
though less common i n models of DSA, m ay also sub- 
stantially heat the ions (IDrurv fc FaHelll986h or electrons 
dGh avamian et al.ll2007l ) of the thermal component. 

A further quantity important for DSA is the efficiency 
of injection of particles from the thermal population into 
the CR population. The fraction of swept up thermal 
protons injected into the acceleration p rocess in SNRs is 
thought to lie in the range 10- 4 -10- 2 (jVolk et alJl2003t 
lEllison fc Cassam-Chena'j2005l) if particle acceleration is 
efficient. 

If the structure of a non-radiative shock is modified 
by CRs, several subtle signature s in the optical emis- 
sion from the shock are expected ([Raymo nd 2001). One 
signature would be a FWHM of the narrow Ha com- 
ponent broader than 20kms~ 1 , the value expected if 
there is no CR acceleration for an upstream medium 
at a temperature of To ~ 10 4 K. Spectra of many 
Balmer-dominated filaments associated with shocks over 
a wide range of Mach numbers, show a narrow Ha com- 
ponent with a FWH M in the range 30-50 km s _1 (see 
ISollerman et al.ll2003l Table 1), indicative of some com- 
mon form of preshock heating. Significant bulk acceler- 
ation of the upstream flow through the precursor would 
also be detectable as a Doppler shift of the narrow com- 
ponent centroid with respect to Ha emission from the 
upstream gas. With the exception of the filament ob- 
served by L07, which also concerns this work, this has 
not been observed yet. Currently, a CR precursor is 
the favored mechan i sm for the inferred preshock heating 
jHester et al.l 11994 ISmith et all 11994 ISollerman et~al] 
120031 : L07), but no self-consistent model of such a precur- 
sor has been compared with data. Here we show that the 
precursor structures predicted by two- fluid models of CR 
modified shocks including particle injection at the sub- 
shock and energy transfer due to the acoustic instability 
are consistent with the above features of Ha spectra. 

Recently, L07 obtained high-resolution Ha echelle 
spectra of an optical filament at knot g in Tycho's SNR, 
covering the postshock region and the ionization pre- 
cursor far upstream. Knot g, located in the eastern 
rim (a = 00 h 25 m 56.5 s , 5 = 64°09'28", J2000.0), is the 
brightest region in Ha emission in the remnant. The 
synchrotro n X-ray emission is also particularly bright i n 
this region (|Decourchelle et aT]|2001l : lHwang et alf 2002). 
Radio data and Hi absorption studies suggest that the 
northeastern rim is decelerating into an inhomogeneous 
amb ient medium, po ssibly the edges of a molecular cloud 
(see lLee et al.| [2004. and references therein). The obser- 
vations by L07 have spatially resolved a steep Ha (narrow 
component) flux increase ahead of the shock discontinu- 
ity, distinct from the photoionization precursor (GOO), 
which L07 attribute to enhanced emission from a CR 
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precursor. They also reported a broadening of the nar- 
row component linewidth by 15kms _1 and a red ward 
Dopplcr shift of the narrow component centroid with re- 
spect to that of the distant upstream Ha emission in this 
region of 5 km s~ 1 . 

In this paper, we provide a self-consistent CR modified 
shock model applied to the observational data from L07. 
We model the shock structure with a time-dependent hy- 
drodynamic two-fluid code. We adjust model parameters 
to obtain a best fit for the calculated spatial Ha profile 
to the observed profile. The two-fluid equations along 
with the free parameters and boundary conditions of the 
shock model are described in Sect. [2 The method of cal- 
culation for the Ha emissivity is given in Sect.[3J A time 
dependent transient solution that provides the best fit to 
the observed spatial Ha profile is presented in Sect.[U In 
Sect. O we discuss the solution space of steady CR mod- 
ified shocks, and propose that the branch of solutions for 
which the CR acceleration efficiency is low, may explain 
the line broadening of the narrow component of the Ha 
line, observed in many SNRs. We discuss our results in 
Sect. El and conclude the paper in Sect. 

2. TWO-FLUID THEORY 

CR modified shocks are sometimes st udied with the 
two-fl uid d escription, firs t deve loped by iDrurv fc V5fkl 
(|1981h and lAxford et all (|1982l ). The CRs are treated 
as a massless fluid exerting a bulk pressure on the ther- 
mal component. The second velocity moment of the CR 
transport equation provides a conservation equation gov- 
erning the CR pressure, which is the quantity that di- 
rectly affects the dynamics of the thermal component. 
Since we are primarily interested in the backreaction of 
the CRs on the gas, it is not essential to follow the CR 
particle distribution in our calculations. 

Two-fluid models have been widely implemented in 
both plane parallel and spher i cally symmetric geom- 
etry (e.g. rMorfill et all flQSl iBoulares fc Coxl 119881: 



Jones fe Kandll992l;lDuffv et al.lli995tlZirakashvili et all 
19961 iFahr et al.l l2000f ) . The basic the ory has been ex- 



tended to incl ude, e .g, wave dissipation ( Volk ct al.l 1984h 
Wagner et a"T1l2007l) . oblique shocks (IWebb et all 119861: 
Frank et al.l Il995f). particle injection ( Kang fc Jonesl 
19901 IZank et all 119931: iKo et alJ Il997t ). and radiative 
cooling (Wag ner et al.l 12006) . The shock structures ob- 
tained with the two-fluid theory agree well with those ob- 
tained from kinetic theory and with those obtain e d wit h 
Monte-Carlo simulations (jKang fc Jonesl 119951 Il997f ). 
The computational expense of solving the two-fluid sys- 
tem numerically is far less than that for the latter two 
methods. 

The approximation of the two-fluid theory consists in 
the limitation that the adiabatic exponent for the CR 
component, 7c, cannot be determined self-consistently, 
and it is therefore usually assigned the constant value 
of 4/3, appropriate for a relativistic gas. In principle, 
the diffusion coefficient for CRs, k, is governed by the 
spectrum of the scattering wavefield, and is thus a func- 
tion of particle momentum, space, and time. Since two- 
fluid models lack information about the CR particle spec- 
trum, a momentum-averaged effective diffusion coeffi- 
cient is commonly used. Although one may prescribe 
the evolution of the closure parameters, n and 7c, (e.g. 
iMarkiewicz et al.lll990t Pones fc Kandll992tlDuffv et all 



Il994h we choose to keep them constant as a first step to 
model the non-radiative shock at knot g. 

The main consequence of the backreaction by CRs on 
the gas is adiabatic heating and compression in the vicin- 
ity of the shock discontinuity where the CR pressure gra- 
dient is largest. When a substantial fraction of the total 
energy has gone into CRs, the overall compression ratio 
exceeds 4 and approaches 7, and the shock discontinuity 
may be entirely smoothed out, i.e. the gas subshock dis- 
appears. Such solutions are referred to as "efficient" or 
"CR dominated" . Conversely, if only a small fraction of 
the shock energy goes into CRs, the solution is termed 
"inefficient" , and the modification of the shock structure 
is weak. 

In the case of fast shocks, such as those responsible 
for Balmer-dominated filaments in SNRs, some regions 
of parameter space permit up to three distinct solutions 
for the same distant upstream conditions and shock pa- 
rameters. Some of these solutions do not exist as time- 
asymptotic steady states, but if exactly two solutions 
exist, one of them is an efficient solution and the other 
is an inefficient one. It is important to perform time- 
dependent runs to determine whether a steady solution 
exists as a time- asymptotic state. 

We neglect the dynamics of the wavefield responsi- 
ble for the scattering of the particles, and thus ignore 
heating of the gas through Alfven wave damping. In- 
stead, we invoke source terms that represent the de- 
cay of sound waves gene rated by the acoustic instabil- 
ity ()Drurv fc Falldll986h . We include injection of CR 
particles at the gas subshock. 

In the remaining parts of this section, we first write 
down the equations for the two-fluid system and the rel- 
evant source terms. We then list the free parameters and 
boundary conditions for our model. 

2.1. Equations 

The following equations govern the two-fluid medium 
consisting of the thermal component and the CR com- 
ponent in plane-parallel symmetry: 



dp dpu _ Q 
dt dx 



dpu dpu 2 
dt dx 



dPn dP, 



c 



dx 



dx 



0, 



(i) 

(2) 



d f pu 2 
di\2~ 



Pa 



1G - 1 



d_ 

dx 



pu 3 JgPqu 
2 7G -1 

-Pop- 
dx 



Pcu 



S G , (3) 



dP c 
dt 



dP c u 
dx 



he 



du 

l)Pc Tx 



d 2 Pc q 



(4) 



x is the spatial coordinate and t is the temporal coor- 
dinate, p, Pa, and T denote the mass density, pres- 
sure and temperature of the gas, and Pc and k are 
the CR pressure and the diffusion coefficient, respec- 
tively. The two fluids move with bulk velocity u. 7c 
and 7c are the adiabatic indices for the gas and of the 
CRs, and are set to the constant values 70 = 5/3 and 
7c = 4/3. Throughout the paper, subscripts G and C 
refer to the gas and the CRs, and subscripts 0, 1, and 2 
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denote a distant upstream value, an immediate presub- 
shock value, and a postshock value, respectively. Equa- 
tion ([T]) expresses mass conservation of the thermal com- 
ponent, and equation ([2]) expresses momentum conserva- 
tion of the thermal component (pu). Equation ([3]) gov- 
erns the total energy density of the thermal component 
{{\/2)pu 2 + Pg/(jg — an d implies that the total en- 
ergy density of the thermal component is conserved in the 
absence of CRs. Equation ^ is derived by multiplying 
the transport equation appropriate for the isotropic part 
of the CR distribution function by (47r/3) p 3 v, where v is 
the CR particle velocity, and integrating that result over 
all CR particle momenta. 

The source terms in equations ([3]) and Q include en- 
ergy transfer from the CR component to the thermal 
component due to the acoustic instability in the CR pre- 
cursor (Sea), as well as CR injection (Sa)' 



Sc — Sci + Sea , 
Sc 



Sr, = — 



ic 



1 



(5) 
(6) 



Both contributions essentially transfer energy from one 
component to the other and the total energy of the two- 
fluid system is always conserved. 

The source terms that represent the energy transfer 
du e to the acoust i c inst ability are identical to those used 
bv lWagner et"afl (|2007l ): 



7c 



1 



5, 



Ca 



Pg[ — 

7G - 1 V a G 
if ■ 



dx 







OPr 



IcPc 



- 1 



dx 



otherwise. 



IcPc 



1 > 0; 



Sound waves are amplified in regions where the rela- 
tive gradient of the CR pressure, jcPc /|dPc/dx|, ex- 
ceeds a critical length scale cig/k- Here, ac = yj JgPg/ P 
is the thermal sound speed. This condition for instabil- 
ity is usually satisfied in the CR precursor because the 
increase in CR pressure scales as k/v s , and v s 3> ag. 
The damping of the sound waves leads to a net energy 
transfer from the CR component to the thermal compo- 
nent. By invoking the source terms (eq.EJ, we bypass the 
initial perturbations and assume that wave damping oc- 
curs at the necessary rate to produce the desired energy 
transfer. The source term is nonzero if the condition for 
acoustic instability is satisfied, r is a time constant that 
determines the rate of energy transfer, and the term in 
brackets ensures that the energy transfer drives the flow 
towards stability. 

We may estimate an order of magnitude value for 
r from dimensional arguments. If we assume a 
Kolmogorov-type turbulence spectrum with the largest 
spatial scale set by the width of the CR precursor, k/v s , 
then the cascade time-scale is approximately 



v s va 



300 yi 



10 24 cm^- 1 
v. 



r)-'(- 



VA 



TOOOkms-M VlOkms-M ' ® 

where va is the Alfven speed. This time scale may be 
comparable to r for CR modified shocks in SNRs. 



I Jones fc Kand (I1990T) fJK 90 hereafter), IZank et all 
(|1993f) . and IKo et alJ l |l997t ) have each adopted a dif- 
ferent approach to simulate particle injection in two- 
fluid models. We follow the method employed by 
JK90, which is a two-fl uid version of that employed by 
iFalle fc Ciddlmral (fl98l : 



Sci — 2 e />l U l (^°G2) 2 W[X - ^subshock) 



0) 



The source term represents the injection of CR parti- 
cles with speed \a,G2 m to the immediate postshock flow, 
and the corresponding rem oval of energy from the ther- 
mal component. Following IFalle fc Giddingsl <| 1987T ) and 
JK90, we set A = 2. The rate at which energy is injected 
is proportional to the mass flux through the shock and 
the parameter e determines the strength of injection. 

Injection is zero ahead of the subshock and has a Gaus- 
sian dependence on distance, w, behind the subshock, 
which is located at ^subshock- The width of the injection 
zone must be at least 15 cells. It must also be much 
narrower than k/v s in order for the solutions not to be 
sensitive to the width of the injection region. Follow- 
ing JK90, injection is turned off smoothly but rapidly 
when the presubshock Mach number with respect to 

ac+Gi = y/haPGi +lcPci)/pi is only slightly larger 
than unity. 

Equations ([TJ> — (J4j) are solved with a second-order finite- 
difference Godunov scheme on a uniform cell grid. The 
CR equation ||3J) is solved implicitly using a Crank- 
Nicholson scheme for the diffusion term. The time-step 
is determined by the Courant-Friedrichs condition with 
the effective sound speed cic+g- The code and technique 
used are identical to those employed by I Wagner et al.1 
((2001 . 

2.2. Free parameters 

Let (f> = Pc I Pg be the ratio of the CR pressure to the 
thermal gas pressure. The free parameters of this model 
are the shock speed, v s (v s — uq in the shock frame); 
the ratio of distant upstream CR pressure to distant up- 
stream thermal gas pressure, 4>q; the distant upstream 
number density, no, and temperature, To; the CR dif- 
fusion coefficient, k; the injection parameter, e; and the 
time scale of energy transfer due to the acoustic insta- 
bility, r. While e, r, and k are essentially unconstrained 
free parameters, we may fix or narrow down the ranges 
of the other parameters from results of previous observa- 
tions. All the above named parameters remain constant 
in time, except for v s which changes as the shock struc- 
ture evolves. 

2.3. Boundary and initial conditions 

For the results presented in this paper, we have as- 
sumed a distance to Tycho's SNR of 2.3 kpc. The non- 
radiative shock at knot g is propagating into a pho- 
toionization precursor, which is at a temperature of 
T re 1.2 x 10 4 K, and n ps 1 cm" 3 (GOO). The 
shock speed inferred from previous shock modeling is 
(~ 20 00±200) kms" 1 (jChevalier et alJlMI iSmith et all 
Il991t G01), although these models have been based on 
the assumption that the available shock energy is con- 
verted solely into thermal energy. The inferred shock 
speed would be higher if a substantial fraction of the 
shock energy went into the CR component. 
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In most environments of the interstellar medium, the 
CR pressure is comparable to the thermal gas pressure 
(|Ferrierd [l998). We have looked at cases for which 
1/3 < 4*0 < 3. In general, we find that the consequences 
of increasing or decreasing (f>o on the structure and evo- 
lution of a shock were very similar to those of increasing 
or decreasing e. We therefore set </> = L 

In the case of knot g, there is some uncertainty in the 
interpretation of the Ha narrow component linewidth, 
due to the presence of an intermediate width component 
(GOO; L07). GOO, and L07 have found that an adequate 
fit to the Ha line profile requires three Gaussian com- 
ponents. The component of intermediate width may be 
produced by protons undergoing secondary charge ex- 
change or by non-thermal motions. Furthermore, the 
assum ption of Gaussian line profiles may not be appro- 
priate ([Raymond et al.ll2008[ ). The presubshock temper- 
ature, Ti, up to which the flow in the CR precursor is 
heated, is therefore taken to be an unknown quantity 
prior to modeling. Our shock model for knot g provides 
an independent estimate of T\ . 

L07 estimated the net acceleration of the flow across 
the precursor, Au, to be in the range 60-130 kms^ 1 . We 
require this condition to be met in our shock model. 

The initial condition (t = 0) for a time-dependent run 
is an ordinary gas shock, for which v s — 2000 km s -1 , 
and for which the CR pressure is homogeneous across the 
grid at the distant upstream value of </>o = 1. The time 
to reach a desired shock structure in any time-dependent 
run should not exceed the age of the remnant, isNR = 
436y. 

We employ free flowing boundary conditions on cither 
side of the grid. 

3. Ha PROFILE CALCULATIONS 

We assume that the optical radiation from the shock 
is not coupled to the hydrodynamics of the flow and cal- 
culate the emission from the system as a separate step, 
using the results from the hydrodynamic simulations as 
boundary conditions. The assumption of separating the 
radiative processes from the hydrodynamics in the CR 
precursor is good because the radiative cooling time is of 
the order of 10 10 s, while the dynamical time-scale across 
the precursor is of the order of 10 8 s. 

In a conventional Balmer-dominated shock, preshock 
neutral hydrogen atoms swept up by the shock (hence- 
forth referred to as neutrals) are initially unaffected 
by the collisionless processes that mediate the shock. 
Within a short distance behind the shock, the neutrals 
may be collisionally excited before being ionized, giving 
rise to a narrow Ha line whose width represents the tem- 
perature of the "cold" preshock neutrals. Preshock pro- 
tons, on the other hand, are instantaneously heated to 
a postshock temperature given by the Rankinc-Hugoniot 
jump conditions. Charge transfer between the hot post- 
shock protons and cold neutrals entering the shock gives 
rise to a population of hot postshock neutrals that are 
responsible for the broad component of the Ha spectral 
line profile. A significant fraction of Ha emission comes 
from Ly/3 trapping; Ly/3 photons from the postshock re- 
gion are converted to Ha photons through scattering. 

However, this conventional picture of a Balmer- 
dominated shock does not explicitly account for a CR 
precursor. The existence of a CR precursor can affect the 



emitted Ha line in a few ways. On the one hand, neu- 
trals are ionized in the precursor and the number of neu- 
trals reaching the shock front is reduced. On the other 
hand, collisional excitation of neutrals within the precur- 
sor serves as an additional source of Ha (and Ly/3) pho- 
tons. The varying temperature and bulk velocity (and 
therefore density) of the precursor flow will also affect 
the rate of Ly/3 trapping. 

We have developed an emission model incorporating 
the above effects of a CR precursor (Lee et al. 2008, 
in preparation), which is briefly described in the follow- 
ing. We consider a flow consisting of hydrogen and he- 
lium. Given a temperature and a velocity profile for the 
precursor, we calculate the spatial evolution of the ion- 
ization states of hydrogen and helium in the precursor 
and postshock regions. The calculations are essentially 
identical to those of G01, except that we have explicitly 
included the effects of the precursor with a given temper- 
ature and velocity profiles. We include collisional ioniza- 
tion by electrons and protons. We distinguish between 
fast and slow neutral hydrogen components, and include 
charge exchange between neutrals and protons. From 
the resulting spatial profiles of slow hydrogen atoms and 
fast hydrogen atoms, we calculate the emissivity of the 
narrow and broad Ha components. For this, we again in- 
clude excitation by both electrons and protons. Since a 
significant fraction of Ha emission comes from Ly/3 trap- 
ping, we similarly calculate spatial profiles of the Ly/3 
emissivity and model the radiative transfer of the Ly/3 
line with a Monte Carlo simulation. The flux of Ha pho- 
tons arising from the radiative transfer of Ly/3 photons 
is added to the intrinsic Ha flux to give the total Ha 
emissivity profile. 

We apply our emission model to the hydrodynamic 
model described in Sect. [5] and compare the results with 
the observed Ha emissivity profile of L07. For this com- 
parison, the distant upstream H neutral fraction is taken 
to be 0.85 and He is fully neutral (GOO). We assume 
equal temperatures for the electrons, ions, and neutrals 
throughout the shock precursor and adopt an equilibra- 
tion fraction between ion and electron temperatures in 
the postshock flow of 0.05. We note that results are 
not very sensitive to the postshock equilibration fraction 
since excitation by protons in the postshock region is 
significant for the assumed shock velocity. 

4. COMPARISON OF MODEL CALCULATIONS WITH DATA 

We have performed a thorough investigation of pa- 
rameter space for the time-dependent two-fluid system 
described in Sect. [2] We find that the observed Ha emis- 
sivity requires the flow in the CR precursor to be heated 
to a presubshock temperature of T\ = 10 5 K. However, 
steady CR modified shock solutions in the high Mach 
number regime are of two types. The solutions tend to 
be either smooth, CR dominated solutions or inefficient 
solutions that are only weakly modified by CRs. In nei- 
ther branch of solutions are the shock structures very 
sensitive to the parameters v s , 4>o, k, r, and e. As a con- 
sequence we could not find a steady solution in which the 
gas in the precursor reached a presubshock temperature 
of Ti = 10 5 K and was accelerated by Am » 1 00 km s" 1 . 
iBlasi et all (|2007f ) and ICaprioli et all (|2008a[ ) explored 
the acceleration time-scale and the evolution of the CR 
acceleration efficiency for SNRs in the Sedov phase, and 
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they found that the shocks evolve towards a quasi-steady 
state on a time scale of the order of 10 3 y. We have, thus, 
found a transient solution that satisfies the observational 
constraints. In the following we first present the shock 
structure of the transient state and then show its subse- 
quent evolution into a CR dominated shock. 



4.1. Transient state solution 

The transient state occurs during the evolution of a 
shock that is initially not modified by CRs to a shock 
that is CR dominated. Figure Q] displays the shock struc- 
ture of the transient solution for which the calculated 
Ha emissivity profiles match the observed profiles in the 
precursor region. In all six panels, the shock profiles are 
shown in the frame comoving with the shock, and the 
flow enters the grid from the right. 

The location of the gas subshock at x = is clear 
from the sudden increase in the broad component flux. 
We take this to be the outermost non-radiative shock 
at knot g, or the outermost part of a corrugated shock 
front (wavy sheet). We are primarily interested in mod- 
eling the Ha narrow component flux upstream of the gas 
subshock. We assume that this emission originates from 
the upstream gas heated by a CR precursor. The model 
must also reproduce the immediate postshock narrow 
and broad component fluxes. The extended emission fur- 
ther downstream is probably due to multiple shock fronts 
superimposed in the line of sight. We do not require our 
model to reproduce the extended emission downstream 
of the immediate postshock region. 

This transient state is reached within 420 y of the 
evolution of a shock for which re — 2 x 10 24 cm 2 s~ 1 , 
r = 426 y, e = 4.2 x 1CT 3 . At time t = 0, the shock 
is not modified by CRs, v s = 2000 kms -1 , and tj> = 1 
throughout the grid. At time t = 420 y (Fig.QJ the shock 
structure is one in which the thermal pressure still dom- 
inates and about 10% of the shock energy has gone into 
CRs. The bottom two panels show that the narrow and 
broad Ha emissivity profiles match the observed profiles. 

The very close similarity between the values of r and 
^SNR is coincidental, although it does reassure us that 
the value for r is physically reasonable. The value of r is 
also consistent with the estimate of the time scale from 
dimensional arguments given in Sect. [2~T1 

The parameters re, r, and e affect the degree of agree- 
ment of the calculated Ha emissivity profile with the 
observed Ha profile in different ways. The spatial ex- 
tent of the observed Ha emission fixes the value of re at 
2 x 10 24 cm 2 s _1 . For a given value of re, we find that r 
primarily governs the balance between heating and de- 
celeration of the gas in the shock precursor during the 
evolution of the shock, while e primarily determines the 
rate at which the CR acceleration (and shock modifica- 
tion) takes place. 

Figure [2] shows the temperature and velocity struc- 
ture of the CR precursor for several shocks that do not 
include injection (e = 0). The model for which the 
values of all parameters other than e are identical to 
those of the model shown in Fig. [T] does not evolve into 
a CR dominated shock. The presubshock temperature 
of the low CR acceleration efficiency steady solution is 
Tj = 3.4 x 10 4 K and Au < 1 kms^ 1 . The calculated Ha 
narrow component profile of this model does not repro- 
duce the observed profile. We find that a lower value of 
t results neither in a larger value of T\ nor a better fit to 
the data. In the absence of injection, the lowest value of 
t for which the solution evolves into a CR dominated so- 
lution is r s=s 698 y. The preshock heating for this shock 
as it evolves towards a CR dominated state, is larger than 
in the case of shocks for which r > 698 y. However, at 
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Fig. 1. — Shock structure of the transient solution for which the calculated Ha emissivity profiles match the observed profiles. The gas 
subshock is located at x = 0. The transient state is reached at time t = 420 y during the evolution of a shock which was initially (t = 0) not 
modified by CRs (see Fig. [3}. The shock parameters are, v a = 2000 km s _1 , cj>o = 1, re = 2 X 10 24 cm 2 s -1 , t = 426 y, and e = 4.2 X 10~ 3 . 
The shock profiles are shown in the shock frame. The upper part of the panel displaying the source functions is the energy transfer due 
to the acoustic instability. The lower part of that panel shows the energy transfer due to injection. The model Ha profile is slightly offset 
along x from the observed Ha profile for visual clarity. 



t — 420 y, this model does not describe the observed Ha 
either. An alternative to invoking injection in order to 
decrease the acceleration time-scale is to chose a smaller 
diffusion coefficient. To obtain a shock structure at 
t = 420 y for which T x = 10 5 K and Au = 100 kms" 1 re- 
quires v s = 2200 kms -1 , k = 10 23 cm 2 s _1 andr = 412 y. 
The calculated Ha profile for this shock model also fails 
to explain the observed profile. The model Ha profiles in 
Fig. [2] are very similar, due to the coarse binning and the 
fact that the emission is dominated by photons from Ly/3 
trapping. The disagreements between the data and the 
results of the models without injection lead us to con- 
clude that injection is a necessary ingredient in a shock 
model that adequately describes the emission from knot 
g- 

The shock structure shown in Fig. [1] can, in fact, be 
reached at an earlier time by increasing the value of the 



injection parameter, e. r need only be somewhat re- 
duced, and k remains unchanged. For example, the evo- 
lution of a shock that is initially not modified by CRs 
at t = will pass through a very similar transient state 
to that shown in Fig. ED at t = 220 y if e = 8.0 x 10~ 3 , 
r = 422 y, and k = 2 x 10 24 cm 2 s" 1 . The condition that 
the time at which the transient state is reached must be 
less than igNR implies that e « 4.2 x 10~ 3 is a lower limit 
for the injection parameter. The transient state cannot 
be reached with a model in which the only CRs are those 
swept up by the shock (i.e. e = 0), even if </>o = 3. 

4.2. Time evolution of model 
In Fig. [3] we show the entire time evolution of a shock 



that acquires the transient state presented in Sect. 14. fl at 
time t = 420 y. The fractional gain in CR pressure ini- 
tially proceeds approximately linearly in time, and the 
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Fig. 2. — Precursor temperature, velocity, and observed and 
model Ha profiles of shocks that do not include injection. The 
different values for v s , t and re are indicated. The normalisation for 
the Ha emissivity profiles is arbitrary and separate for the narrow 
and broad component profiles. The model profiles are slightly offset 
along x with respect to each other and the observed profile for 
clarity. The model Ha profiles are very similar, due to the coarse 
binning and the fact that the emission is dominated by photons 
from Ly/3 trapping. None of the shock models without injection 
give a satisfactory fit to the Ha emissivity data (particularly the 
narrow component). We conclude that injection is necessary in a 
model that can explain the data. 

transient state at t — 420 y occurs during a phase of 
evolution in which the shock modification is still rela- 
tively weak. The transition from a weakly CR mod- 
ified shock that contains a gas subshock to a smooth 
CR dominated shock is quite rapid. The shock speed 
recedes to v s w 1830 km s -1 . During this phase a den- 
sity "spike" develops due to a temporary overcompres- 
sion arising from the combined compression across the 
CR precursor and the gas subshock. Within a region of 
width comparable to the diffusion length-scale, the den- 
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Fig. 3. — Evolution of a shock, which was initially not mod- 
ified by CRs, and for which v 3 = 2000 km s , <f>o = 1, 



2 X 10 24 cm 2 s~\ r = 426 y, e = 4.2 X 10 -3 . The transient state 
shown in Fig. [T] occurs during a phase in which the thermal pres- 
sure dominates. During the rapid transition to a CR dominated, 
smooth shock, beginning at t ~ f 000 y, a density spike forms (see 
Fig. 2). All profiles are shown in the frame comoving with the 
shock at t = 0. The time-space diagram shows the same results as 
the sequence of profiles. The flow enters the grid from the top. The 
vertical line in the time-space diagram and in the panel showing 
the fraction of shock energy going into CRs as a function of time 
marks the instant in time the transient state occurs (t = 420 y). 

sity is enhanced by a factor greater than 20 with respect 
to the distant upstream density, and a factor of 3 with 
respect to the density in the immediate postshock region. 
The density spike travels downstream at approximately 
210 kms -1 with respect to the smooth CR dominated 
shock. In Fig. [4] we show the density, velocity, temper- 
ature and pressure structure of the shock after it has 
developed into a CR dominated shock (t = 1600 y in the 
evolution shown in Fig. [3]). 

The density spike in tim e-dependent CR modified 
shocks was first noticed bv iDorfiJ (11984D . and its for- 
mation was explained bv iDrur-yT i 1987f ) and JK90. The 
portion of the postshock flow in which the CR pres- 
sure dominates is bounded by the forward shock and the 
density spike. The flow upstream of the density spike 
is supersonic with respect to oq but subsonic with re- 



Cosmic-ray precursor of a Balmer-dominated shock 



9 



Flow density and speed 



2000 



1000 i 




Fig. 4. — Shock structure at time t = 1600 y in the evolution 
of the model shown in Fig. \3\ A density spike has formed during 
the rapid transition to a CR dominated shock and is now traveling 
downstream. The flow upstream of the density spike is supersonic 
with respect to ag, but subsonic with respect to ac+G> while the 
flow downstream of the spike is subsonic with respect to og. The 
density spike may be long-lived and may play a role in the enhance- 
ment of the magnetic field downstream and affect the morphology 
of synchrotron emission from shocks that accelerate CRs. 

spect to ac+Gi an d the flow downstream of the den- 
sity spike is subsonic with respect to ac- The den- 
sity spike may t herefore be long-live d as it moves down- 
stream, although lJun fe Jones! (|1997l ) have demonstrated 
that the flow through the density spike forms Rayleigh- 
Taylor instabilities that reduce the density enhancements 
slightly. Its formation certainly warrants attention in 



spherically symmetric or 3D simulations of SNR blast 
waves as it may play a role in the enhancement of mag- 
netic fields and the morpho logy of synchrotron emission 
(|Cassam-Chenai et al"1l2007l ). 

A factor of 2-4 enhancement in the thermal X-ray 
brightness with respect to the local average may be a 
signature for a density spike travelling downstream of a 
strongly CR modified shock. This feature would proba- 
bly be difficult to detect due to projection effects, con- 
tamination from ejecta clumping near the forward shock, 
and low temperatures of the thermal component (see 
Fig.SJ. 

5. WEAKLY CR-MODIFIED STEADY SHOCKS AND THE 
Ho LINEWIDTH OF BALMER-DOMINATED FILAMENTS 

The solution space of two-fluid models for steady, 
high Mach- number (M > 30), adiabatic, CR modified 
shocks contains solutions of high CR acceleration effi- 
ciency (i.e. (f>2 > 1) and low CR acceleration efficiency 
(i.e. 4>2 <C f). For some combination of shock parame- 
ters and di stant upstream conditions there exist multiple 
solutions. iBecker fe Kazanai (|2001l ) derived exact ana- 
lytic expressions for the domains in solution space that 
contain multiple solutions, and the domains in wh i ch the 
solutions exhibit a gas subshock. iMalkov et al.l (|20Q0T ) 
have investigated the conditions under which such a bi- 
furcated system may self-regulate. 

If Sg = Sc = and if (j>o is less than some critical 
value (f> c that depends on the shock Mach number, then 
for any given distant upstream state there are three pos- 
sible postshock states. Of the three solutions, the one for 
which th e CR acceleration effi ciency is intermediate, is 
unstable (|Mond fe Drurvill998l). and does n ot exist as a 
time- asymptotic state ( Donohue et al.l fl994l ) . The solu- 
tion for which CR acceleration is most efficient exists if 
(and only if) it is smooth, i.e. it does not contain a gas 
subshock. 

Finite source terms determine whether a shock that 
is accelerating CRs evolves into a CR dominated shock, 
or a shock for which <p2 remains small. For high Mach- 
number shocks 4> c <C 1 if Sg — Sc = 0. In this case, 
the solutions bifurcate only if the distant upstream CR 
pressure is low compared to the pressure of the thermal 
component, and steady solutions are mostly CR dom- 
inated shocks. However, if Sea < 0, we find that (f> c 
approaches unity, i.e. the solution space is bifurcated if 
the distant upstream CR pressure is comparable to the 
thermal gas pressure. This gives rise to the existence of 
the low CR acceleration efficiency branch in high Mach- 
number shocks for which </>n ~ <f> c ~ 1- 

In Fig. [5] we show the structure of a shock for which 
v s = 2000 kms" 1 , c = 2 x 10 24 cm 2 s~\ (f> Q = 1, r = 
634 y, and e = 1.0 x 10 -5 . The pair of values for r and 
e lie near the upper boundary of values beyond which 
a low CR acceleration efficiency solution does not exist 
(see Appendix). The presubshock temperature reaches 
T\ = 4.2 x I0 4 K and the net acceleration across the 
precursor is Au = lOkrns -1 . 

The shock structures in the low CR acceleration effi- 
ciency branch of steady solutions do not vary much with 
shock speed in the range 300 kms -1 < v s < 3000 kms -1 , 
which we have explored. We find that in these shocks 
no more than 1% of the shock energy is channelled into 
the CR component. The choice of t = 634 y and e = 
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Fig. 5. — Profiles of a steady shock in the low CR acceleration efficiency branch of solutions. v s = 2000 km s -1 , ft = 2 X 10 24 cm 2 s~ 1 , 
4>o = f, t = 634 y, and e = 1.0 X 10 — 5 . The gas subshock is located at x = 0, and the profiles are shown in the frame of the shock. The 
upper part of the panel displaying the source functions is the energy transfer due to the acoustic instability. The lower part of that panel 
shows the energy transfer due to injection. The set of values for r and e yield the maximum bulk acceleration of the flow through the CR 
precursor, Am m lOkrns -1 . The temperature in the precursor reaches T\ = 4.2 X 10 4 K immediately ahead of the gas subshock. Most 
solutions in the low CR acceleration efficiency branch are shocks structures similar to the above. 

1.0 x 10 -5 (Fig. [5]) yield the largest value for Au that is 
possible for a steady shock in the low CR acceleration ef- 
ficiency branch of solutions, for which v s = 2000 km s^ 1 , 
k = 2 x 10 24 cm 2 s _1 , and 4>o = 1- 

For very small values of r and large values of e, we 
find solutions in which the presubshock temperature ex- 
ceeds 10 5 K and Au — > 0. In the greater part of param- 
eter space, however, including cases for which e = 0, the 
preshock temperature reaches 2-6 x 10 4 K, and Au = - 
lOkrns -1 . We also note that the steady solutions of the 
low CR acceleration efficiency branch are reached very 
quickly, usually within t < 100 y, if the initial conditions 
at t = are those of a shock that is not modified by CRs. 

The insensitive nature of the low CR acceleration effi- 
ciency branch of steady solutions to r, e, and v s , pro- 
vided that they do not exceed critical values for low 
CR acceleration efficiency steady solutions to exist, may 
be an explanation for the small range of FWHM (SO- 
SO km s -1 ) for the narrow component of the Ha line 
observed in Balmer-do minated filaments of many SNRs 
(|Sollerman et al.|[2003T ). Currently, other models of non- 
radiative shocks which do not include CRs are not able 
to predict the FWHM of 30-50 kms -1 . Furthermore, the 
lack of observed bulk Doppler shift in the narrow com- 
ponent is also consistent with the small values for Au 
that we obtain. Since only a small fraction of the shock 
energy goes into the CR component, the shock speeds 
inferred from previous models applied to non-radiative 
shocks in SNRs remain valid. 

6. DISCUSSION 



The results from the time-dependent solutions pre- 
sented in Sect. S] depend on the assumed distance to Ty- 
cho's SNR. The canon i cal va lue, also assumed here, is 
d = 2.3 kpc. iLee et all ()2004f ) have briefly reviewed the 
debate on the distance to Tycho's SNR. While most au- 
thors adopt d = 2.3 kpc as established bv lChevalier et all 
( 1980) and subsequen t ly confirmed by severa l other stud- 
ies (see [Strom 19881) . iSchwarz et al.l (|1995l ) argue for a 
value of 4.6 ± 0. 5 kpc, based on 21 c m ab sorption fea- 
tures. However, iBlack &: Raymond! (|1984D found that 
absorption features at those velocities were present in 
the spectra of stars within 2.5 kpc tha t were close to Ty - 
cho on the sky. The calculations by IVolk et all (2007) 
together with the high energy 7-ray flux upper limit 
from HEGRA {Ah aronian et alj feOOl) imply a distance 
4 kpc > d > 3.3 kpc. If the distance to the remnant is 
greater than 2.3 kpc, then the shock speed inferred from 
proper motions, and the intrinsic width of the CR pre- 
cursor would be larger. The values for n and r must be 
scaled accordingly to obtain a model whose calculated 
spatial Ha emissivity profile matches the observed pro- 
file, e must also be increased to allow for the desired 
transient state to be reach ed within t < isNR- However, 
Ivan Adelsberg et all (|2008l ) have applied their improved 
shock models to the Ha data and obtained a 15% slower 
shock speed of v s — 1600-1700 kms -1 under the assump- 
tion that Pc/Pq = 0. The slower shock speed would 
lead to a correspondingly smaller distance to Tycho's 
SNR. Our result that ~ 10% of the shock energy goes 
into CRs during the current state of the shock at knot g 
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would imply a 5% increase in v s , pa rtially compensating 
for the lower shock speed inferred by lvan A dclsbcr g et al.l 
(l2008h . 

Although the spatial emissivity profile of the Ha nar- 
row component places a stringent constraint on k, a sepa- 
rate upper limit for the width of the CR precursor comes 
from the fact that a substantial fraction of hydrogen 
atoms must avoid ionization ahead of the gas subshock in 
order to produce the observed broad component. Thus, it 
is important that our model describes the narrow compo- 
nent flux from the CR precursor as well as the immediate 
postshock broad to narrow component flux ratio. For the 
values of the diffusion coefficient, distant upstream den- 
sity, and required preshock temperatures in this problem, 
the CR precursor length-scale is much shorter than the 
preshock electron, proton, and photon ionization length- 
scales, and the H ionization fraction remains nearly con- 
stant throughout the CR precursor. This justifies our 
choice to adopt the distant upstream neutral fraction of 
85% derived from the photoionization model by GOO as a 
fixed boundary condition. A different choice of upstream 
ionization fraction would affect the immediate postshock 
ratio of broad to narrow component fluxes, as well as 
the rate of Ly/3 trapping in the CR precursor and, there- 
fore, the spatial Ha profile in front of the gas subshock. 
The upstream density no only plays a small role in de- 
termining the Ha narrow component flux from the CR 
precursor through enhancing Ly/3 trapping in the precur- 
sor region close to the gas subshock. The more significant 
property that determines the spatial flux profiles is the 
shock temperature profile, in particular the value of the 
presubshock temperature. 

The presubshock temperature T± = 10° K inferred 
from our model is larger than the temperature of ~ 
4 x 10 4 K implied by the narrow component linewidth 
in the data of L07. A possible reason for this is the 
single-temperature approximation we use for the precur- 
sor gas. The neutrals may have a lower temperature if 
the length scale for charge exchange is a significant frac- 
tion of the CR precursor length-scale. This may be the 
case for shocks propagating at speeds v s > 2000 kms -1 
into a medium of low fractional ionization. 

It should be noted that the calculated Ha emissivity is 
sensitive to the ratio of electron to ion temperatures in 
the CR precursor. We have adopted a ratio of unity, with 
the assumption that some plasma-physical process pro- 
vides the requisite electron heating throughout the CR 
precursor. Possible processes include the resonant ex- 
change of energy between electrons and protons via lower 
hybrid waves excited, for example, by the two-stre am in- 
stability due to shock-reflected ions (|Lamingll200u) . or by 
the same mecha nism that induces high frequency mag- 
netosonic waves (|Ghavamian et al.ll2007l ). 

In the mechan i sm fo r wave dissipation proposed by 
iGhavamian et all |2007) the electrons may attain a tem- 
perature of up to 0.3 keV in a CR precursor. For a shock 
for which v s = 2000 km s -1 and the preshock density is 
uq = Iran -3 , this would imply a precursor width less 
than ~ 10 16 cm to avoid complete ionization of the neu- 
trals. This is smaller than the precursor width observed 
in knot g, suggesting that the electron temperature is 
considerably lower. 

Based on a study of the Tycho re mnant's X-ray 
morphology and spectral characteristics, IWarren et al.l 



(2005) proposed that the forward shock at most az- 
imuthal angles is strongly CR modified. A shock com- 
pression ratio approaching 7 could explain the proximity 
between contact discontinuity and forward shock, though 
projection effects may allow a smaller compression ratio 
(jCassam-Chenai et al.l 12008). Knot g is situated at an 
azimuthal angle approximately 80° east of north, and 
is clearly recognizable as a local minimum in the ratio 
of the radius of the cont act discontinuity to that of the 
forward shock (Fig. 4 of I Warren et al.|[2005f) . Assuming 
that the X-ray shock and the Balmer-dominated shock 
at knot g are directly associated, the small ratio of radii 
in the direction of knot g is consistent with our shock 
model in the transient state: the shock in the transient 
state is still in the phase in which the shock modification 
due to CRs is relatively weak, although it is developing 
into a CR dominated shock. 

Although the evolution of the shock in our time- 
dependent run into a CR dominated state (Sect. 14. 2p 
is physical within the framework of the two-fluid the- 
ory (Sect. (2), the question of whether it is a realistic 
description of the subsequent evolution of the Balmer- 
dominated shock at knot g must be approached with 
care. Our plane-parallel model accounts neither for adia- 
batic losses due to expansions perpendicular to the shock 
normal, nor for particle escape upstream. These two 
effects may reduce the CR acceleration efficiency ap- 
preciably. FOT_e2camj3k, the ste ady state, kinetic shock 
models of lCaprioli et al.l (|2008al ) that include particle es- 
cape, and the spheri cally symmetric, kinetic models of 
iKang fc Jones! (|2006f ) predict that the acceleration effi- 
ciency, which they define as PC2I Po u Oi reaches 60% dur- 
ing the Sedov phase, whereas we obtain an equivalent effi- 
ciency greater than 70% in our time-dependent run. The 
absence of losses in our model should not significantly 
affect the early phase of the evolution during which CR 
modification is still weak, and during which the transient 
state that describes the Balmer-dominated shock at knot 
g occurs. 

The shock models presented in this work are based on 
the interpretation of the observations by L07 that the 
steep rise in the flux of the Ha narrow component ahead 
of the shock front is due to a CR precursor. The pos- 
sibility of other types of precursors remains. Of these, 
a fast neutral precursor is currently thought to be the 
most likely alternative candidate. A fast neutral pre- 
cursor is mediated by the hot postshock neutrals (those 
responsible for the broad component of the Ha line), 
which escape upstream and deposit some energy via 
ch arge exchange and elast ic collisions. T he calculations 
bv lLim fc Ragal ()1995l ) and lKorreckl (|2005l ). however, pre- 
dict that the net heating due to a fast neutral precursor is 
too small to account for the observed narrow component 
line broadening. In contrast to a CR precursor, the effi- 
ciency of preshock heating by fast neutrals is sensitive to 
the degree of thermal equilibration between electrons and 
ions, the shock speed, and the neutral fraction upstream. 
It is therefore difficult to explain the narrow range in line 
broadening seen in many Balmer-dominated filaments in 
which diverse shock conditions obtain. 

The possibility that the narrow component flux in- 
crease ahead of the gas subshock is due to a superpo- 
sition of multiple shocks in the line of sight has not been 
ruled out. However, as L07 have argued, a superposi- 



12 



Wagner et al. 



tion of shocks would imply a gradual flux increase of the 
broad component ahead of the gas subshock in the same 
manner as the flux of narrow component. This is not 
observed. 

7. CONCLUSIONS 

In summary, CR acceleration in the forward shocks 
of SNRs results in the heating and acceleration of the 
preshock medium which may explain some features of 
the optical emission of Balmer dominated filaments. We 
have found a transient state in the evolution of a shock 
from one that is initially not modified by CRs to one 
that is CR dominated, for which the calculated Ha 
emissivity profile matches the emissivity profile across 
the Balm e r-dom inated filament in knot g observed by 
iLee et aD (|2007l ). The values of the parameters for this 
shock model are an initial shock speed v s = 2000 km s , 
a distant upstream CR pressure to thermal gas pressure 
ratio (f)Q = 1, a diffusion coefficient k — 2 x 10 24 cm 2 s _1 , 
an energy transfer time-scale due to the acoustic instabil- 
ity r = 426 y, and a lower limit to the injection parameter 
e = 4.2xl0~ 3 . 

The structure of steady shocks that belong to the low 
CR acceleration efficiency branch of solutions for fast 



shocks are relatively insensitive to the values of r, e, 
K, and v s in the range 300 kins" 1 < v s < 3000 kins -1 , 
provided that the parameters are chosen such that a 
steady solution in the low CR acceleration efficiency 
branch exists. The solutions are usually reached as time- 
asymptotic states within less than 100 y, even if e = 0. 
The mild heating of the preshock gas up to typically 2- 
6 x 10 4 K, and the negligible bulk acceleration of the flow 
in the precursor (Au < 10 kms" 1 ) may provide a natural 
explanation for the characteristic broadening of the nar- 
row component linewidth that is observed to lie in the 
small range FWHM = 30-50 kms -1 in many SNRs, and 
the lack of bulk Doppler shift of the narrow component 
observed for these cases. 
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APPENDIX 

In the following we give an example of the limits for r and e within which the low acceleration efficiency branch of 
steady CR modified shocks, as described in Sect.[5j exist. 

If e = 0, the exact value for <fi c depends only on the dimensionless quantity r\ = tv^/k. Th e method to find the 
range in r\ for which two distinct solutions exist for a given distant upstream state is described in lWagner et al.l ()2007f ) 
for radiative shocks of lower Mach number. Here, we find, for example, that for shocks for which v s = 2000 km s -1 , 
<f)Q = 1, K = 2 X 10 24 cm 2 s -1 , and e = 0, the low CR acceleration efficiency branch of solutions exists up to r a rj 634 y, 
and the (smooth) high CR acceleration efficiency branch of solutions exists down to 77, ~ 63 y. In the range tj, < t < r a 
both an inefficient and an efficient steady solution are allowed. 

As expected, a large injection parameter drives the solution towards one that is CR dominated. For example, if 
v s = 2000 kms -1 , k = 2 x 10 24 cm 2 s _1 , <fio — 1, and r = 426 y, the critical value for e above which a low CR 
acceleration efficiency solution cannot exist is ~ 1.0 x 10~ 4 . Conversely, a shock for which v s , k, <po are the same as 
above but e = 4.2 x 10~ 3 (the value obtained from the model in Sect. [4]) can only remain in the inefficient branch of 
solutions if r < 63. 
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